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Forced Convection in High 
Porosity Metal Foams 

This paper reports an experimental and numerical study of forced convection in high 
porosity (e~~0.89-0.97) metal foams. Experiments have been conducted with aluminum 
metal foams in a variety of porosities and pore densities using air as the fluid medium. 
Nusselt number data has been obtained as a function of the pore Reynolds number. In the 
numerical study, a semi-empirical volume -averaged form of the governing equations is 
used. The velocity profile is obtained by adapting an exact solution to the momentum 
equation. The energy transport is modeled without invoking the assumption of local 
thermal equilibrium. Models for the thermal dispersion conductivity, k d , and the inter- 
stitial heat transfer coefficient, h sft are postulated based on physical arguments. The 
empirical constants in these models are determined by matching the numerical results 
with the experimental data obtained in this study as well as those in the open literature. 
Excellent agreement is achieved in the entire range of the parameters studied, indicating 
that the proposed treatment is sufficient to model forced convection in metal foams for 
most practical applications. [S0022- 148 1(00)01903-4] 
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Introduction 

Forced and buoyancy-induced convection in porous media have 
been studied extensively for over 50 years (see Kaviany [1] for a 
good review on the subject). However, most studies in this cat- 
egory have been restricted to packed beds and granular materials, 
since they have direct application to naturally occurring porous 
media, with porosities in the range 0.4-0.6. As discussed below, 
there are relatively few investigations of transport phenomena in 
very high porosity media (e>0.9) such as metal foams. 

Only during the last ten years, transport phenomena in metal 
foams have started to receive attention ([2-6]). Hunt and Tien [2] 
studied forced convection in metal foams with water as the fluid 
phase. Using the technique of volume averaging, and under the 
simplifying assumption of local thermal equilibrium, they showed 
that a single energy equation could adequately describe forced 
convection in metal foams. Sathe et al. [3] studied combustion in 
metal foams as applied to porous radiant burners. Lee et al. [4] 
studied the application of metal foams as high-performance air- 
cooled heat sinks in electronics packaging. In their experimental 
study, they demonstrated that a 1 cm 2 chip dissipating 100 W 
could be cooled using an aluminum foam heat sink, and a low- 
power muffin fan. Recently, Calmidi and Mahajan [5] addressed 
heat conduction in aluminum metal foams. They obtained experi- 
mental data and developed an analytical model for the thermal 
conductivity, taking into account the open-celled structure of 
metal foams. Bhattacharya et al. [6] extended that work to cover a 
wider range of solid to fluid conductivity ratios. Further, the ana- 
lytical mode! was refined to address issues of isotropicity. 

Our review suggests that there are no thermal transport studies 
in the literature (pertaining to high porosity metal foams) that (a) 
employ experimentally validated thermal conductivity data, (b) 
consider local thermal nonequilibrium among phases, and (c) ac- 
curately consider metal foam microstructure. The focus of this 
paper is to take these factors into account. 

Metal foams (Fig. 1) can be broadly classified as porous media 
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in which the medium has a distinct but continuous and rigid solid 
phase, and a fluid phase. They are typically available in high 
porosities (e>0.9). They also have a unique open-celled structure. 
As a result, most of the past studies on packed beds, and granular 
porous media are not applicable to metal foams. There is thus a 
need to develop a fundamental understanding of the transport phe- 
nomena in these foams. As shown later in this paper, both numeri- 
cal and experimental studies are currently needed to generate heat 
transfer information for practical applications. 

Energy transport in porous media has generally been studied 
with the assumption of local thermal equilibrium among the two 
phases. That is, a single homogeneous equation is used to describe 
transport with the assumption that the solid and fluid phases have 
the same temperature field. Following such an assumption, the 
effective stagnant thermal conductivity, k e , of a porous medium 
is used to account for the solid and fluid phase conductivities. The 
stagnant conductivity k f , is sometimes augmented by the disper- 
sion conductivity, k d . Dispersion is a hydrodynamic phenomenon 
which becomes prominent at high Reynolds numbers, especially if 
k e is small in magnitude. Several empirical and analytical studies 
have attempted to quantify dispersion in porous and fibrous media 
([2,7-9]). The general conclusion is that at high Peclet numbers, 
k d is linearly proportional to the local flow velocity. 

Under certain situations, however, the assumption of local ther- 
mal equilibrium has to be relaxed. Such a condition might arise 
when there is internal heat generation ([3]) or when the difference 
in the stagnant conductivities of the solid and the fluid phases is 
significant. The latter situation is encountered when metal foams 
are used with low-conductivity fluids like air or water. With local 
thermal nonequilibrium, the solid and fluid temperature fields 
have to be obtained separately. In addition, appropriate coupling 
is necessary between the solid and fluid phase energy equations to 
account for the interfacial heat transfer. 

The focus of this study is forced convection in metal foams. 
First, the experimental method used to obtain Nusselt number data 
for a number of aluminum metal foam samples with air as the 
fluid medium is described. Following a discussion of the appro- 
priate form of the governing equations of momentum and energy 
in metal foams, the details on the numerical method used to solve 
them are provided. Finally, the numerical results are compared 
with the experimental data to obtain values of empirical constants 
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Fig. 1 Picture of a typical metal foam sample 



in the thermal dispersion and the interfacial heat transfer coeffi- 
cient models. 

Experiment 

A schematic of the experimental setup is shown in Fig. 2. Alu- 
minum metal (alloy T-6201) foam samples of size 63 mmX45 
mmX196 mm were used in the experiments. Table 1 lists the 
characteristics of the samples that were used in the experiment. 
The porosity is the void volume fraction and the pore density is 
the number of pores per unit length of the material (PPI, pores per 
inch) and, d f and d p are the average fiber diameter and pore 
diameter, respectively. The porosity of each sample was estimated 
using the weight of a given volume of the sample and the density 
of aluminum alloy T-6201. The fiber diameter, d f , was measured 
using a microscope, and the pore diameter, d p , was estimated by 
counting the number of pores in a given length of material. Both 
d f and d p are average values. The pore density (PPI) is a nominal 
value supplied by the manufacturer. 

For each experimental run, a metal foam sample was placed in 
a plexiglass tube. Air flow was achieved by connecting the plexi- 
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Fig. 2 Schematic of the experimental setup used for the 
forced convection experiments 



glass tube to. a fan/motor assembly (not shown in figure) down- 
stream of the test sample. The direction of air flow is from left to 
right as shown in Fig. 2. Very low conductivity styrofoam was 
used to block the space around the sample to fully prevent flow 
bypass. A pressure transducer was used to measure the pressure 
drop across the metal foam sample. An orifice plate was used to 
measure the flow rate through the metal foam sample. 

The metal foam samples were brazed to 9.5-mm thick alumi- 
num skins on two sides. On one of the skins, grooves were cut to 
isolate a section of the metal foam. Patch heaters were fixed to 
this isolated portion of the skin between the grooves (Fig. 3). The 
resulting physical dimensions of the heated section of the metal 
foam are 63 mmX45 mmX 1 14 mm. The grooves in the skin were 
loosely filled with insulating material to prevent conduction losses 
from the heated to the unheated sections. Then, two thin (2-mm) 
grooves were cut in the metal foam below the grooves on the skin 
(Fig. 3). This was done to prevent back-conduction of heat from 
the heated to the unheated sections. This isolating scheme ensures 
that fully developed flow enters the heated section of the metal 
foam (Calmidi [10]). Six 36-gage T-type thermocouples (axial 
thermocouples) were fixed along axial locations in the skin as 
shown in Fig. 3. For some samples, two thermocouples were fixed 
at off-center locations (25 mm from the center). Two additional 
thermocouples were used to monitor the ambient temperature and 
the temperature of the air at the midplane before it entered the 
heated section. The total flow rate was estimated according to 
ASME/ISO-5167-l:1991 using the pressure drop across the ori- 
fice plate as well as the temperature and pressure of the air enter- 
ing the orifice plate. The average velocity (superficial) of air flow- 
ing through the metal foam, u 0 , is then the average flow rate 
divided by the area of cross section. 

During a typical experiment run, the speed of the motor was 
adjusted such that the flow velocity was set at a desired value. The 
power input to the patch heaters was set between 15-35 W. The 
temperature readings indicated by the thermocouples were moni- 
tored until they reached steady state. This took approximately 
5-10 minutes depending on the flow velocity. It typically took 
longer at lower flow velocities. The temperatures were monitored 
for an additional interval of five minutes to ensure that steady 
state had indeed been reached. Steady state was assumed if the 
temperature did not vary more than 0.1 °C during a five-minute 
interval. The temperature readings of all the thermocouples with 
respect to the ambient temperature were noted. The temperature 
variation along the axis (^-direction) was within 15 percent of the 
average value, while the temperature variation perpendicular to 
the axis ^-direction) was less than 0.1 °C in all the experimental 
runs. The air temperature before it entered the heated section 
(T in ) was also measured. It was slightly higher than the ambient 
temperature (7^), but always less than two percent of the maxi- 
mum wall temperature, indicating that back-conduction effects 
were negligible. Figure 4 shows the data (wall temperature varia- 
tion) collected from a typical experiment for Sample 2 (Table I): 
The data collected for all the samples is documented in Calmidi 
[10]. 

Using the measured wall temperature data, the heat flux is con- 



Table 1 Characteristics of metal foam samples used in experimental study 
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The effect of varying the power input to the heaters was also 
investigated to verify that nonlinear effects like buoyancy, radia- 
tion, and property variation were negligible. In all cases studied, 
the Nusselt number was found to be nearly independent of the 
input power. 

The experimental error estimate can be made based on the ac- 
curacies in the measurement of the individual quantities. The error 
(thermocouple calibration and resolution of the data acquisition 
device) in the estimation of A7\, vg is 0.3°C Based on the pub- 
lished accuracies of the voltmeter and ammeter, the error in the 
estimate of q is two percent. The errors in the estimate of the 
length and area are neglected because they are extremely low 
(<0.3 percent of the measured quantity). Based on these values, 
the error in the estimate of the Nusselt numbers is 4.3 percent. 
This estimate is based on a representative value of A7" avg of 8°C. 
Based on the accuracy of the pressure transducers, the error in the 
estimate of the velocity is 4. 1 percent. 



Fig. 3 Schematic of side and top views of metal foam samples 
used in the experiments 



verted into an average Nusselt number. The average temperature 
Ar avg is defined as the average of the six values of the tempera- 
ture measured. That is, 




n - T h 



(1) 



where ai is the number of points where the wall temperature was 
measured. It is noted that, although this value is not the actual 
average wall temperature, it does represent an appropriate scale to 
represent the temperature difference. Next an average heat trans- 
fer coefficient is defined as 

h = q/{A±T avg ) (2) 

where q is the power input to the heaters and A is the area of cross 
section of the metal foam. From this, an average Nusselt number 
is calculated using 



Nu = hL/k e 



(3) 



where L is the length of the heated section of the metal foam, and 
k e is the effective thermal conductivity (k/ e + k 5e ) as given in 
Table 1. 



Analysis 

The steady volume-averaged momentum equation that governs 
fluid flow in porous media is given by 



— y V -uu ~ — V p + -V 2 u- — u- -==||u||u. 



(4) 



In Eq. (4), the drag due to the fibers is represented by the familiar 
non-Darcy formulation. The third term on the R.H.S. of Eq. (4) is 
the traditional Darcy term which is used to account for the pres- 
sure drop due to viscous friction at the walls of the fibers, and the 
fourth term is the nonlinear correction and represents the pressure 
drop due to the form drag and the flow separation that takes place 
at higher Reynolds numbers. A linear superposition of the two 
terms as shown here is semi -empirical, but has been known to 
match experimental data in many situations ([2,8,11-13]). A 
modified form of Eq. (4) was used by Calmidi [10], Chapter 3, to 
determine the permeability, K, and the inertia coefficient,/, of the 
metal foam samples. These are listed in Table L 

The energy equations for the solid and the fluid phases based on 
thermal nonequilibrium are 
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Fig. 4 Variation in wall temperature with average flow velocity (Sample 2, 
Table 1) 
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+ h 5f a sf (T-T f ). (6) 

The coupling between the two equations is achieved using the 
interfacial term which represents the heat transfer between the two 
phases via a heat transfer coefficient, h sf , and the specific surface 
area, a sf . It is important to note that Eqs. (5), and (6) are semi- 
empirical and are not the full volume-averaged forms of the en- 
ergy equations ([14,15]). The full volume-averaged treatment is 
considerably complex, and has been solved exactly only for a 
one-dimensional system of parallel capillary tubes ([16]). How- 
ever, it is noted that the semi-empirical treatment employed here 
has been successfully used in the past to model transport in beds 
of spheres ([17,18]). 

An analytical model for obtaining the effective conductivity of 
metal foams has been derived in ([5]). Their experimentally vali- 
dated model was based on the structure of metal foams. It is 
reproduced as 
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where r = 0.09, and the ratio bIL is given by 
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The solid-phase conductivity, k se , is obtained by setting k f - 0, in 
Eq. (7). By doing so 
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Similarly, k fe is obtained by setting £,=0 in Eq. (7), 
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Values of and £ /f obtained using Eqs. (9) and (10) are listed in 
Table 1. 

In Eqs. (5) and (6), the dispersion conductivity is assumed to be 
isotropic. That is 



Based on earlier work ([2,7]), the following model is proposed for 
the dispersion conductivity: 



k d u 
-^(Re^Pr,)-. 



(12) 



In Eq. (12), the coefficient of thermal dispersion, C D , remains to 
be determined. As a simplification, this model for k d does not 
account for wall effects ([19,20]) except through the change in the 
velocity profile near the wall. 

One of the most comprehensive models for the interfacial heat 
transfer coefficient, h sfy for packed beds is by Wakao et al. 
([17]). For foamed materials, however, no such general model 
exists. However, the radial temperature gradients within the fibers 
are expected to be small and it is reasonable to use the Nusselt 
number correlation for external flow over a body of an appropriate 
cross section. Based on a correlation developed by Zukauskas [21] 
for cylinders in crossflow in the range of Reynolds numbers (40- 
1000), a model for the interfacial heat transfer is proposed as 
follows: 



Nu s/ =^ = C r Re 05 Pr a37 . 



(13) 



In Eq. (13), an empiricism has been introduced in the form of the 
undetermined coefficient C 7 , The Reynolds number is based on 
the cylinder diameter and Pr is the fluid Prandtl number. In apply- 
ing this correlation, the noncircular cross section of fibers has to 
be considered. However, upon examination of the Nusselt number 
correlations for external flow over bodies of varying cross sec- 
tions ([22]), it appears that the exponent of the Reynolds number 
has a weak dependence on the cross section. The multiplicative 
constant in Eq. (13), however, varies considerably. Based on these 
observations, Zukauskas* correlation has been modified: as shown 
in Eq. (13). Further, it is noted that, in Eq. (13), the Reynolds 
number is based on the fluid velocity near the fiber. Hence, due to 
the presence of the solid matrix, its point value is Kt-{ud f iev). 
The solid-fluid interfacial surface area for arrays of parallel cylin- 
ders intersecting in the three mutually perpendicular directions is 



lird f 



(14) 



For metal foams, this expression is modified by taking the struc- 
ture into account (i.e., open cells shaped like dodecahedra, and 
noncircular fiber cross-section). By doing so, the values of d p and 
d f , as given in Table 1, are multiplied by 0.59, and 1 
-e~ ((l ~ e)/0 04) , respectively (see [10], Chapter 1) before being 
used in Eq. (14). 

The transport equations (Eqs. (4)-(6)) are nondimensionalized 
using -X=x/H, Y=yIH, U=ulu 0 , 0=(T~ r /n )/A7 avg where 
u 0 is the average velocity entering the channel, H is the height of 
the channel, A7 avg is the average wall temperature, and T in in the 
inlet fluid temperature. The resulting equations are 
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(18) 

Note that Bi, and are not constants but vary along the 
^-direction, since h s j (Eq. (13)) depends on the velocity. In Eqs. 
(15)-(18), all quantities except C D and C T are known. 

The heat dissipation from the heated wall occurs both through 
the fluid phase and the solid phase. Hence, the total heat trans- 
ferred can be written as 



avg I 
JO 



El 



h sf a $f Hy[K a specified error of 10~ 5 . This method of calculating the velocity 
£ . profile was checked by solving the full two-dimensional momen- 
tum equation, Eq. (4) in a rectangular domain with an aspect ratio 
(length to height) of two using the ADI Method of Peaceman and 
Rachford (see [24]). The fully developed velocity profile obtained 
from this simulation was identical to the exact solution given by 
Eq. (21). Since computational time could be saved considerably 
by using the exact solution, it was employed for all further 
computations. 

The energy equations were solved using the ADI finite differ- 
ence scheme of Peaceman and Rachford ([24]). The convective 
term in the fluid-phase equation was discretized using first-order 
upwind differencing for increased stability. The diffusion terms in 
the solid-phase and fluid-phase equations were discretized using 
three-point stencils in each direction. To avoid energy balance 
problems, the numerical domain was extended by five percent at 
_ hi _ f UH lk se d$ s k fe dd f \ the inlet and the exit. Extensions of 10 percent and 15 percent had 

Nu ~ ~ e ~~ \ VJT W + 7" 17 1' ^ 20) no effect on * e lem P erature fie l d or the Nusselt number in the 
J 0 e e ' heated section. The equations were marched to asymptotic steady 

state. The convergence criterion for steady state was that the 
change in the solid-phase and fluid-phase temperatures was less 
than 10" 4 between successive iterations. At steady state, a heat 
balance within one percent was obtained. A time-step of size 0.03 
was found to be an efficient value. While reducing the time-step 
below this value did not have any effect on the final solution, the 
convergence rate slowed down. 



Writing Eq. (19) in nondimensional form, an average Nusselt 
number for the heated surface is 



Equation (20) has been formulated to directly compare experi- 
mental results with model predictions. 

Numerical Procedure 

A schematic of the numerical domain is shown in Fig. 5 along 
with the orientation of the coordinate axes. Numerical boundary 
conditions are also shown in Fig. 5. For the heated wall, the ex- 
perimentally determined wall temperature profile is specified. The 
temperature gradients for both the solid phase and the fluid phase 
are zero on the insulated wall. At the inlet, the fluid is assumed to 
enter at the inlet temperature T in . At the exit, a zero diffusion 
condition is assumed for the fluid phase. This assumption is rea- 
sonable, since advection is expected to be the dominant mode of 
transport in the streamwise (x) direction. The solid phase is as- 
sumed to be insulated at both the inlet and the exit. 

A fully developed velocity profile is prescribed at the inlet. This 
velocity profile (function of Y alone) is calculated using the exact 
solution for the second-order nonlinear ODE, Eq. (15). Following 
Vafai and Kim [23], the exact solution due to flow through a 
channel is given by expressing the pressure drop in terms of the 
center-line velocity. It is reproduced below: 



1-1 — — l(cosh(D(l-2K+C))) 2 



4 = e/Re/6Da 2 £ = e/(4Da 2 ) + 2A 
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(21) 



Since the average velocity u 0 , and not the center-line velocity is 
known, the center-line velocity had to be guessed and iterated 
until mass balance was achieved. That is, the integral of the non- 
dimensional velocity over the height of the channel is unity within 



Grid Dependency and Code Validation 

A constant grid spacing was used in the X-direction and a vari- 
able grid spacing was used in the K-direction. The variable grid 
was laid out in three layers. In Layers I and III, adjacent to the 
heated and insulated walls, respectively, a constant grid spacing 
was used. In Layer II, sandwiched between Layers I and III, a 
variable grid was used. Table 2 summarizes the results obtained 
during the grid dependency tests for one of the test samples (see 
[10], Chapter 4) with Re*=84.1. 

The first column in Table 2 refers to the constant grid spacing 
A Y in Layers I and III in the K-direction. The second and third 
columns refer to the number of grid points in Layer II and the 
total number of grid points in the K-direction, respectively. The 
fourth column is the number of grid points in the heated section 
along the X-direction. Based on the grid dependence study, the 
grid parameters listed in the second row were used for all 
computations. 

In order to validate the code, numerical results were compared 
to those of Yokoyama and Mahajan [25] as summarized in Table 
3. The numerical results ([25]) are indicated within parentheses. 
Pe* in Table 3 is defined according to ([25]). 



Table 2 Results from grid dependency simulations 
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Fig. 5 Schematic of numerical domain with boundary condi- 
tions 
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Table 3 Comparison of numerical 
Yokoyama and Mahajan [25] 


results with those in 
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Results and Discussion 

Using the numerical code, computations were performed for the 
conditions of the experiments to verify if the results could be 
matched in a reasonable manner, and to determine appropriate 
values for C D and C r . The experimental data was collected in 
the range of nondimensional parameters, Re^l 1-135, Da 
~5X10"" 3 -1X10~ 2 , and Pr e ~2.6x 10' 3 -7.5x 10 -3 . Note 
that the effective Prandtl number values are extremely small com- 
pared to the fluid Prandtl number (0.707) due to the high thermal 
conductivity of the metal foam. 

First, the condition of zero dispersion (i.e., C D = 0) was studied. 
Figure 6 shows a plot of the Nusselt number as a function of 
Re^Pr, for different values of C T with C D = 0. The results are for 
Sample 2, and the experimental points are shown by the solid 
symbols. As expected, the Nusselt number increases with increase 
in C T since the efficiency of heat transfer between the two phases 
improves. However, note that there is practically no difference in 
the computed results for C r >0.52 indicating that the phases are 
close to the thermal equilibrium condition. Figure 7 shows the 
variation of the temperature profile along the transverse direction 
at the midplane for three different Reynolds numbers. Note that 
the adiabatic wall temperature (i.e., the wall temperature at Y 
= 1 ) for the solid and fluid phases is not the same as the ambient 
temperature for all three values of Re*. However,- as Re* in- 
creases, the adiabatic wall temperature approaches the ambient 
temperature. The indication is that the thermal boundary layer at 
the heated wall decreases with increase in Re* . 



To understand the relative magnitudes of the various transport 
phenomena, it is useful to examine the solid and fluid-phase heat 
transfer rates at the heated wall. In Fig. 8, the local variations of 
the nondimensional heat transfer rates (the two components of Eq. 
(20)) are plotted. Clearly, the fluid-phase component accounts for 
a negligibly small portion of the total heat rate. Further, its local 
value decreases rapidly along the streamwise direction. This can 
be explained as follows. The predominant mode of transport from 
the heated surface to the fluid is by conduction through the solid 
phase and interfacial heat transfer from the solid to the fluid 
phase. Due to the efficient interfacial heat transfer, the fluid tem- 
perature is close to the solid temperature. This has an effect of 
decreasing the fluid temperature gradient at the wall (by thicken- 
ing the boundary layer), and consequently the local heat transfer 
rate. Further, since the fluid conductivity is much smaller than the 
solid conductivity, the net heat transferred from the wall to the 
fluid is low. 

In the numerical results presented so far, thermal dispersion has 
not been taken into account. It can be seen from the plot in Fig. 8 
that the heat transferred directly to the fluid from the heated wall 
is negligible compared to the heat transferred interfacially from 
the solid to the fluid. Hence, the dispersion phenomenon that pri- 
marily enhances the wall heat transfer, is expected to be low. 
Figure 9 shows a plot of the total Nusselt number as a function of 
Re^Pr, for Sample 2. The plots are for two different values of C T 
as shown in the figure. For each C T value, two lines are shown. 
They are the results of numerical simulations with and without the 
inclusion of dispersion effects. For the simulations for which dis- 
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Fig. 8 Solid and fluid phase local heat transfer rate. Sample 2, 
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Table 4 
Tien [2] 



Characteristics of foam samples used In Hunt and 



No. 


Material 


e 


d p (cm) 


A'MO 7 
(m 2 ) 


/ 


(W/m-K) 


(W/m-K) 


HI 


C 


0.97 


02 


4.1 


0.11 


0.073 


0.58 


H2 


C 


0.97 


0.1 


0.92 


0.077 


0.073 


0.58 


H3 


Ni 


0.97 


0.1 


0.96 


0.089 


0.854 


0.579 


H4 


Al 


0.97 


0.2 


4.8 


0.17 


2.19 


0.579 


H5 


Al 


0.94 


0.25 


17.0 


0.30 


3.84 


0.558 




Fig. 10 Nusselt number as a function of Re^r, for different 
values of C D (Sample H3) 



persion effects were included, a value C D = 0.1 was used. As it is 
evident from the figure, dispersion has no significant effect on the 
total transport. This is immediately clear by re-examining the ther- 
mal dispersion model postulated earlier, 



^-C 0 (Rc r Fr # ). 



(22) 



For k d to have any significant effect on the overall transport, it 
must be at least of the same order of magnitude as k e . For a 
typical value of Re^Pr f = 0.3, and C D = 0.1, k d is merely three 
percent of the total stagnant thermal conductivity. Although its 
value is still seven times higher than the fluid conductivity, the 



effect on the overall transport is negligible because, as mentioned 
before, the dominant mode of transport is by conduction through 
the solid phase and interfacial convecuon, which is largely unaf- 
fected by dispersion. The value, C o = 0.1, used for illustrating this 
is not entirely arbitrary. It is based on order of magnitude esti- 
mates in ([7]). In any case, it is clear that it is difficult to accu- 
rately quantify dispersion effects when the solid-phase conductiv- 
ity is much larger than its fluid-phase counterpart, as is the case 
here. 

Although an accurate estimate of C D cannot be made by using 
the metal foam-air experiments, an estimate of C T is possible. 
From Fig. 9, it is seen that a good match of the data with the 
numerical results is obtained with C r =0.52. This value is repre- 
sentative of heat transfer from cylinders in crossflow ([21]). Fur- 
ther, as will be seen later, this value of C T is sufficient to model 
forced convecuon transport in metal foams with water as the fluid 
phase. We note here that a two-equation model for energy may 
not be necessary for all combinations of experimental data. How- 
ever, a detailed study that establishes bounds for nondimensional 
parameters for which two-equation analyses are necessary may be 
useful for application to practical situations (see, for example, 
[26]). 

The only other experimental and theoretical work that addresses 
forced convecuon in metal foams is that of Hunt and Tien [2]. 
Their study used foams made of aluminum, nickel, carbon, with 
water as the fluid phase. Hence, the fluid conductivity is not a 
negligible fraction of the effective conducuvity. Consequendy, 
dispersion effects are expected to be much higher. Although an 
error was found in the numerical calculations (which eventually 
led to the conclusion that C D - 0.025), the experiments seemed to 
be in order. Of the seven samples used in their experiment, five 
representative samples were chosen for analysis here (see Table 
4). The conductivities in the last two columns of Table 4 were 
evaluated using the expressions developed in this study (Eqs. 
(9),(10)). Since water was used as the fluid phase, interfacial heat 
transfer was efficient and all samples satisfied local thermal equi- 
librium conditions as expected. (C 7-= 0.52 was used for the com- 
putations, based on the aluminum-air study described earlier.) 

To study the dispersion effect numerically, Sample H3 was 
chosen for investigation. The average heat transfer coefficients 
were extracted from their experimental Nusselt numbers and used 
here for comparison. The computed Nusselt numbers for different 
values of the dispersion coefficient, C D , are shown in Fig. 10 
along with the experimental data. A reasonable fit is obtained for 
C D = 0.06. As expected, and unlike in the aluminum foam-air 
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Fig. 15 Nusselt number as a function of Re* for the 40 PPI 
samples (Table 1) 
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combination used in this study, the dispersion effect is consider- 
able here. In this particular case, heat transfer enhancement due to 
dispersion is almost 80 percent. Using Eq. (22), and a representa- 
tive value of Re^Pr e = 30, the dispersion conductivity is calculated 
to be 175 percent of the stagnant conductivity. This is a significant 



enhancement compared to the low value of three percent obtained 
earlier for the aluminum-air combination. The experimental re- 
sults in ([2]) are plotted in Fig. 11 for a thermal dispersion coef- 
ficient value of 0.06. A good fit is obtained for all samples, sug- 
gesting that C D ~ 0.06 is appropriate. 

By analyzing the data collected in this study for foam- air com- 
binations, and that published in Hunt and Tien [2] for foam-water 
combinations, the coefficients C T and C D values of 0.52, and 0.06 
seem reasonable. As a last check, numerical simulations were per- 
formed for the same experimental conditions of all samples (Table 
1) used in the present study. The results of these simulations are 
plotted in Figs. 12-15 as a function of the Reynolds number, 
along with the experimental data. Once again, a very good fit with 
the experimental data is obtained confirming the values of C D and 
C T to be good. However, the following point is noted. There 
appears to be a systematic deviation between the experimental and 
computed results for high Reynolds number. In particular, the 
numerical results appear to underpredict the experimental Nusselt 
number for Re^>80. The indication is that some other transport 
enhancing effect (e.g., turbulence) may be taking effect. No such 
systematic deviation is observed in the experimental results of 
Hunt and Tien [2]. However, the maximum Reynolds number in 
that study was around 60. 

Summary 

A detailed study of forced convection in metal foams has been 
performed. The goal of the study was to quantify thermal disper- 
sion and thermal nonequilibrium effects in metal foams. To this 
end, both experimental and numerical methods were employed. 
Experiments were performed with a wide variety of aluminum 
metal foams. Results indicate that for foam-air combinations, the 
transport enhancing effect of thermal dispersion is extremely low 
due to the relatively high conductivity of the solid matrix. How- 
ever, for foam-water combinations, ours and past available results 
indicate that thermal dispersion can be very high and accounts for 
bulk of the transport. 

A thermal nonequilibrium model was used in all cases. The 
values C T =0.52, and C o = 0.06 were found to be appropriate. 
Currently, forced convection transport in foamed materials is be- 
ing studied using foams of other materials (carbon, copper) along 
with various fluids (air, water, FC-72) in order to reconfirm the 
results of the present study. Those findings will appear in a future 
publication. 
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Nomenclature 

a sf = specific solid-fluid interfaciaJ surface area (nT 1 ) 
B '/ = h 5f a sf H^KIk fty fluid-phase effective Biot number 
Bi, = h SJ a s/ H 2 /k se . solid-phase effective Biot number 
C D = coefficient of thermal dispersion 
C T = coefficient in Eq. (13) 
Da = yjKIH, Darcy number 

/ = inertia coefficient 
h sf = interfaciaJ heat transfer coefficient (W/m 2 -K) 

K = permeability (m 2 ) 

= t*c p /k e , Prandd number based on effective conduc- 
tivity 

q = heaunput to patch heaters (W) 

= u 0 ylK/v. Reynolds number based on permeability 
T in = measured fluid temperature (before entering heated 
section) (°C) 

T'amb = measured fluid temperature (before entering unheated 

section) (°C) 
T wi = measured wall temperatures (°C) 
AT avg = average wall temperature with respect to fluid inlet 
temperature, T in 
u Q = measured average flow velocity entering foam sample 
(m/s) 

Greek Symbols 

e = (r-r in )/(A7" avg ) 

Subscripts 

d = dispersive 
e - effective 
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